Critical behavior of plastic depinning of vortex lattices in two dimensions: Molecular 

dynamics simulations 
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Using molecular dynamics simulations, we report a study of the dynamics of two-dimensional 
vortex lattices driven over a disordered medium. In strong disorder, when topological order is 
lost, we show that the depinning transition is analogous to a second order critical transition: the 
velocity-force response at the onset of motion is continuous and characterized by critical exponents. 
Combining studies at zero and nonzero temperature and using a scaling analysis, two critical expo- 
nents are evaluated. We find u ~ (F - Fcf with /3 = 1.3 ± 0.1 at T = and F > F^, and v ~ T^'^ 
with — 0.75 ±0.1 at _F = Fc, where Fc is the critical driving force at which the lattice goes from 
a pinned state to a sliding one. Both critical exponents and the scaling function are found to exhibit 
universality with regard to the pinning strength and different disorder realizations. Furthermore, 
the dynamics is shown to be chaotic in the whole critical region. 

PACS numbers: 74.25. Uv, 64.60. Ht 
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I. INTRODUCTION 



Due to the competition between interactions and ran- 
domness, the dynamics of coherent structures driven over 
a disordered medium exhibit a great variety of phases. In 
particular, the transition from a pinned state to a slid- 
ing one occurring at a critical driving force Fc^ known as 
the depinning transition, is both a great theoretical chal- 
lenge and relevant for numerous systems: superconduc- 
tor vortices, colloids, Wigner crystals, magnetic bubbles, 
charge density waves, magnetic domain walls... It has 
been suggested on phenomenological grounds by Fisher— 
that the depinning transition could be regarded as a crit- 
ical phenomenon in which the velocity and the driving 
force would respectively be the order parameter and the 
control parameter. Though it was originally intended to 
describe the depinning of charge density waves (CDW), 
this idea has proven very useful in many other domains. 
In the elastic limit, it has been shown theoretically that 
for most manifolds and for CDW the depinning does be- 
have like a second order transition with a power law re- 
sponse V ^ [F — FcY at the onset of motion and /3 < 1—. 
When it comes to situations in which an elastic descrip- 
tion is no longer valid, however, the theoretical descrip- 
tion of the phenomenon is much more difficult. In partic- 
ular, the nature of the depinning transition (continuous 
or discontinuous) remains an open problem. A continu- 
ous depinning transition (second order) with /3 > 1 is ob- 
served in experiments^^— and numerical simulations^^—. 
On the other hand, experiments on CDW—, and 3D nu- 
merical simulations of vorticesi^ suggest a discontinuous 
(first order) depinning transition where the velocity-force 
curve displays hysteresis and jumps between pinned and 
unpinned states. On the theoretical side, various models 
have been developed to describe non elastic dynamics. 
A coarse grained model has been proposed^^ in which a 
visco-elastic coupling is used as an effective description of 
topological defects or phase slips. In the mean field limit. 



it predicts two kinds of depinning: a continuous one, 
belonging to the universality class of elastic depinning, 
and an hysteretic one. The existence of an hysterctic 
depinning in a special case of this model has been con- 
firmed using functional renormalizationM. Other phase- 
slip models also predict hysteresis for CDW—. Con- 
versely, numerical studies of a model focusing on the 
filamentary nature of the flow in 2D show a continuous 
depinning with (3 ~ 1.5^^. Other phase-slip models based 
on the XY model also suggest the absence of hysteresis 
at the thermodynamic limitl^ (in 3D) and a continuous 
second-order plastic depinning transition with an expo- 
nent j3 ~ 1.7— (in ID). Another approach is to introduce 
inertia in the equation of motioi*i^. In this case a con- 
tinuous depinning transition is found for small incrtial 
parameter, otherwise a discontinuous transition with hys- 
teresis is found. In periodic systems with a displacement 
field of dimension N = 2 {e.g. superconductor vortices, 
colloids, Wigner crystals), simulations have established 
that strong disorder leads to dislocations and plasticity. 
At the depinning threshold, regions of pinned particles 
(zero velocity) coexist with particles flowing around the 
pinned region a^i^i^"" — . Besides, the transition seems to 
be continuous and smooth (/3 > 1) in 2D (d = 2 is here 
the dimension of the embedding space)2ri£. Simulations 
on 2D colloids also indicate that the relaxation time near 
depinning obeys a power law, as expected in a second or- 
der transitions^. However, most studies in which critical 
exponents are evaluated are carried out in the case of a 
ID displacement field (N — 1, e.g. Josephson junction 
arrays, metallic dots, CDW), while there arc very few 
such studies for N ~ 2. In particular, there is to the best 
of our knowledge no available study on superconductor 
vortices giving the /? exponent in the 2D plastic regime 
(it should be noted though that /3 has been numerically 
evaluated for superconductor vortices in 3D2i). 

In this paper, we perform molecular dynamics simula- 
tions of 2D vortex lattices with strong random pinning 
and study the depinning transition induced by an exter- 
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nal driving force. Our system belongs to the category 
of 2D periodic systems with a 2D displacement field and 
short range interactions. It could model 3D supercon- 
ductors (either conventional or layered) in an effective 2D 
regime, i.e. when the vortex line tension is high enough 
for the lines to remain straight. 

The behavior of the system near the depinning transi- 
tion is studied at both zero and nonzero temperature, 
which allows an accurate measurement of the critical 
driving force. At zero temperature, the depinning is con- 
tinuous and highly plastic. The static channels observed 
at the onset of motion are identified with the so-called 
single particle regime and considered as a finite size ef- 
fect. In the intermediate range of driving force, the mo- 
tion is chaotic, whereas at high driving force the lattice 
reorders and chaos disappears. In the threshold vicinity, 
a study of the temperature dependence of the velocity 
allows us to determine the true critical force despite the 
single particle regime. Moreover, it is shown that the ve- 
locity scales as a power law of temperature vf=f^ ^ T^^^ 
at F = Fc and the associated critical exponent 5 is de- 
termined. At T = 0, a second power law v ~ [F ~FcY is 
found with an exponent /? > 1. Varying simultaneously 
the driving force and the temperature, we then find ev- 
idence of the existence of a scaling law, which confirms 
the values of /? and 5 evaluated respectively at T = and 
F = Fc. This analysis has been performed for various 
system sizes, disorder realizations and pinning strengths, 
indicating that both /3, S and the scaling function exhibit 
some degree of universality. 

II. NUMERICAL MODEL 

As in Ref. [13, we study Ny Abrikosov vortices interact- 
ing with Np random pins in the {x, y) plane. We consider 
the London limit Ai ^ ^, where Al is the penetration 
length and ^ is the coherence length, i.e. we treat vortices 
as point particles. The overdamped equation of motion 
of a vortex i at position Vi reads 

j^i p (1) 

where rip is the distance between the vortex i and the 
pinning site located at Tp, r^j is the distance between 
the vortices i and j located at and r^, and Vi is the 
2D gradient operator acting on r,;. rj is the viscosity 
coeflScient. = FSq is the Lorentz driving force due 
to an applied current. F**^ is the thermal gaussian white 
noise with zero mean and variance 

{Ft^^it)F;!lit')) = 2^kBT5,j5^J{t - t') 

where fijiy = x,y and ks is the Boltzmann constant. The 
vortex- vortex repulsive interaction is given by a modified 
Bessel function 



and the attractive pinning potential is given by 

where Rp the radius of the pins, and a^, and ap are 
tunable parameters. Depending on the relative strengths 
of the vortex-vortex and vortex-pin interactions, the dy- 
namics can be either dominated by elasticity or disorder. 
The strength of the vortex-vortex interaction is fixed 
by setting ay ~ 2.83 10~'^ALeo where cq is an energy 
per unit length. The relative disorder strength ap/ay 
is then chosen high enough for the depinning transition 
to exhibit plasticity (in our model, plasticity is found 
above ap/ay ^ 0.01). Molecular dynamics simulation is 
used for Ny ~ 30n^ vortices in a rectangular basic cell 
{Lx,Ly) = (5, 6-\/3/2)riAL, with n from 3 to 8, i.e. Ny 
from 270 to 1080. Periodic boundary conditions arc used 
in both directions. The vortex-vortex interaction is dealt 
with using a neighbor list method with a cutoff radius 
Tc ~ 7.1 Al- The number of pins is set to Np = 5Ny, 
and their radius is Rp = 0.22 Al. The average vortex 
distance is oq = Al. We use a unit system in which 
77 = 1, Al = 1, eo = 1 and ks ~ 1. 



III. BEHAVIOR AT r = 

We start with a perfect lattice at high velocity and 
slowly decrease the driving force until the system reaches 
a pinned state. The force is then varied back and forth 
with various force steps in order to check for hystere- 
sis. The whole process is done at T = 0. In Fig. [1] we 
plot the typical shape of the average velocity-force curve, 
showing four distinct regions. These regions are better 
illustrated by the typical trajectories of vortices^^ dis- 
played in Fig. [2l 

Just above the depinning threshold (region I), the flow 
occurs along one or several non crossing static channels 
while the rest of the system remains pinned (see Fig. [2]). 
The velocity-force curve exhibits jumps and hysteresis. 
This is related to the existence of several sets of such 
channels for a given driving force. Each set corresponds 
to an hysteresis branch while jumps in velocity are in 
fact jumps from a branch to another one, i.e. from a set 
of channels to another one. Note that the branch cho- 
sen by the system depends on the force ramping rate. 
Moreover, on a given channel, the positions of vortices 
are fully determined by the position of one of them. In- 
deed, it can be verified by choosing one vortex at a given 
position on the channel, and plotting the positions of all 
the other vortices. Since the channel is one dimensional, 
this means that the dynamics can be modeled by a sin- 
gle degree of freedom. Because of the periodic boundary 
conditions, this single degree of freedom sees a periodic 
potential, leading to the so-called single particle regime: 
the system experiences a saddle-node bifurcation at the 
critical force ; above F^ the velocity is a periodic func- 
tion of time and scales as v ^ {F — FcY^'^. We measure 



3 




Figure 1. (color online) Schematic of the average velocity v 
versus driving force F. The vertical dashed lines separate the 
different kinds of flow observed: pinned static lattice (region 
0), single particle regime (region I), disordered chaos (region 
II), smectic chaos (region III) and decoupled channels (region 
IV). The dotted line crossing the horizontal axis at _F = F* is 
the prolongation of the curve from region II to region I (see 
text section HV A|) . 



indeed this power law on each branch corresponding to 
a given configuration of the channels. Moreover, simula- 
tions of boxes shaped as long strips in the longitudinal 
direction show that the range of force in which region 
I is observed decreases when the longitudinal size is in- 
creased. At large sizes, the size of region I seems to shrink 
to zero. As a result, we assume that region I vanishes in 
the infinite size limit, and should not be taken into con- 
sideration for the study of the critical behavior. 

In region II and III, vortices flow along changing 
interconnected channels (see Fig. [2]). The transition 
from region II to region III is defined by the appearance 
of transverse smectic order indicated by small peaks in 
the structure factor along the ky axis, while the sketch 
of the high velocity channels becomes visible. It occurs 
near the inflexion point in the velocity-force curve, also 
known as the peak in differential resistance in reference 
to the experimental tension-intensity curves. In region 
IV, vortices stop switching channels (see Fig. [2]) and a 
linear behavior v ^ F — Fc is observed. 

Following Ref. [22l . we now show that the motion in 
regions II and III is chaotic by evaluating the sign of 
the maximum Lyapunov exponent A. The existence of 
at least one positive Lyapounov exponent illustrates the 
sensitive dependence on initial conditions, which is a sig- 
nature of chaos. In Fig. [3l we plot the distance 
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in the phase space between two neighboring trajectories 
{v\{t),...,Y\^{t)) and (r2(t),...,r^^(t)). Two distinct be- 
haviors are found. In regions II and III, d{t) grows ex- 
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Figure 2. (color online) Typical trajectories of the vortices 
in the four moving regions shown in Fig. [T] single particle 
regime (region I), disordered chaos (region II), smectic chaos 
(region III) and decoupled channels (region IV). 



ponentially at first and then saturates. The exponential 
growth indicates the existence of at least one positive 
Lyapunov exponent, which proves the existence of chaos. 
A saturation effect appears when d{t) becomes of the or- 
der of the size of the chaotic attractor. On the other 
hand, in regions I and IV d{t) remains constant at large 
times (the largest Lyapunov exponent is zero) indicating 
a closed orbit in phase space while the decrease at short 
times is due to the transient regime associated with the 
existence of negative Lyapunov exponents. 



IV. CRITICAL EXPONENTS AND SCALING 

In the original approach of Fisher—, it is argued that the 
depinning transition is a critical phenomenon implying 
scaling laws and critical exponents near threshold. It is 
therefore expected that vf=f^ ^ T^/"^, and vt=o,f>f^ ~ 
{F- Fc)^, where S and /3 arc critical exponents. In the 
following, we show that our data support the existence of 
such a critical phenomenon: the critical exponents S and 
/? are determined, and a scaling relation is found between 
the velocity, temperature and driving force. Note that 
the threshold behavior that we observe in our simulations 
is effectively continuous since hysteresis and jumps have 
only been observed in region I, which we believe vanishes 
in the infinite size limit. 
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Figure 3. (color online) Distance d between two neighboring 
trajectories versus time in the region I, II, III and IV. 

A. Critical force and exponent 5 

First of all, we want to determine the critical force F^- 
The main issue is that because of hysteresis in region I, 
there is no unique dcpinning force Fc directly accessible 
to measurements. Consequently, we have to evaluate an 
effective critical force F* by prolongating the velocity- 
force curve from region II to region I. The intersection 
with the w = axis defines F* (see Fig. [T|). This evalu- 
ation can be improved by studying the temperature de- 
pendence of the velocity close to the depinning threshold. 
As shown in Fig. 21 two different behaviors arc observed, 
depending on which side of the transition the force is. 
Above the critical force, v approaches a non zero limit as 
T goes to 0, leading to convex curves with an horizon- 
tal asymptote on the left in logarithmic scale. Below the 
critical force, v goes to faster than a power law, result- 
ing in concave curves in logarithmic scale. The change 
in convexity when the force is varied indicates that the 
effective critical force F* has been crossed. The results 
obtained by these two methods are consistent, and com- 
bining both allows to improve the accuracy and precision 

of^^;. 

Furthermore, in agreement with a second order phase 
transition, we can extrapolate aX F = F* a power law 
rcsponssi (a linear response in logarithmic scale at the 
change of convexity) from which we measure the critical 
exponent 5: 

VF=F. ^ T^" (3) 

As we shall sec in section ITV CI this is consistent with the 
existence of a scaling relation, and the value of 5 obtained 
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Figure 4. (color online) Average velocity v versus temperature 
T at low temperature for different forces around the critical 
force F* (from bottom to top: 10^ x F = 11; 12; 12.5; 13; 14), 
OLp/a-u = 1.05 and A''„ = 270. The dashed black line is the ex- 
trapolation of v(T) at the critical force from which we extract 
. Inset: value of 5^^ for different samples: 2 different dis- 
order strengths ap/a^ = 1.05 and a^ja^ = 0.35, 3 different 
sizes (N^ = 270 for ap/a„ = 1.05, = 270, 480 and 1080 for 
cep/ctv — 0.35), and different sets of positions of the pins (2 
sets for ap/a^ — 1.05 and = 270, 5 sets for ap/a^ — 0.35 
and Nv = 270). The error bars correspond to the different 
lines one can draw to extrapolate a power law behavior at the 
change of convexity. 

here allows to collapse all data available on a single curve, 
supporting the validity of extrapolating a power law be- 
havior at = F* . This analysis has been performed on 
9 samples with different pinning strengths {ap/ay = 1.05 
and 0.35), different system sizes (iV„ = 270 to 1080) and 
different realizations of the random positions of the pin- 
ning centers. The resulting values of are shown in 
the inset of Fig. H leading to = 0.75 ± 0.1. 

B. Critical exponent /3 

We now go back to the T = case to study the 
velocity-force curve using the values of F* obtained in 
the previous section. The mean velocity v versus effec- 
tive reduced force f = {F - F*)/F* is plotted in Fig. [5] 
for ap/ay = 0.35, showing the existence of a power law 
regime in the critical region, which lies from the lower 
boundary of region II (/ ~ 0.1) to close to its upper 
boundary (/ ^ 1) 

t'T=o,/>o (4) 

This power law results in a linear region in logarithmic 
scale, which slope (3 has been measured on the 9 samples 
of section IIV Al No significant differences between the 
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Figure 5. (color online) Average velocity v versus effective 
reduced driving force f = {F — F*)/F* (circles and squares), 
including the critical region (squares). Bottom right inset: 
value of P for different samples (same samples as Fig. [H 
Qfp/Q^ = 1.05 and 0.35, iV„ = 270 to 1080 and different sets 
of positions of the pinning centers). Top left inset: values of 
F* for different samples (left side: ap/av = 1.05, 2 samples ; 
right side: OL^ja^ = 0.35, 7 samples). 



samples were measured, leading to a unique value /? = 
1.3 ±0.1 (sec bottom right inset of Fig. [5]). The precision 
of /? is limited by the precision of F* . Note that using 
only the T = results (see section IIV Ap to determine 
the critical force gives a similar result, except that the 
uncertainty on /3 is greater. 



C. Scaling law 

The power law dependence of v versus both / and T 
strongly suggests to go on with the identification of the 
depinning transition with a critical phenomenon and to 
look for evidence of a scaling relation between the veloc- 
ity, driving force and temperature. First of all, we want 
this relation to be expressed in terms of dimensionless 
quantities, and to be independent of the prefactors in 
the two power laws mentioned before. Wc thus define 
and To such as 



W/>0,T=0 = Wo/'' 



T 



1/5 



(5) 



Considering dimensionless velocity v = v/vq and tem- 
perature T = T/Tq, we define the scaling function S as 
follows: 



v\f\ 



s± (r|/| 



-135 



(6) 



where the unknown branches 5*+ and of the scal- 
ing function correspond respectively to / > and 



/ < 0. Moreover, the observed power law dependences 
'^/>o,T=o = E^nd Vf=o = T^/* imply that S±{x) satis- 
fies 



lim S+(x) 



1 



lim x-^/^S±{x) = 1 (7) 



Graphically, it means (in logarithmic scale) that 5+ is 
asymptotic to the horizontal axis for T — )■ (driving 
dominated regime) while both 5+ and S- have an oblique 
asymptote with slope for / — > (temperature domi- 
nated regime). The intersection of these two asymptotes 
occurs at a; = 1, defining a force dependent crossover 
temperature T ~ To\f\^^ between the two regimes. Note 
that the change of variables (w,T) — > (w, T) is equiva- 
lent 
orig 



10' 



10 



vlff 



10 



10' 



10" 




10 



10 

TIfl 



10 



(55 



Figure 6. (color online) Scaling plot «|/|~^ versus Tl/l"^^ 
with /3 = 1.3 and = 0.75 featuring the 9 samples shown 
in the insets of Fig. [4] and Fig. [5] {ap/a^ = 1.05 and 0.35, 
A*'„ = 270 to 1080 and different sets of positions of the pinning 
centers). 

We check for the existence of the scaling by plotting 
in Fig. [6]'!;|/|~^ versus T\f\~^^. We observe a collapsing 
of data to a single curve (with two branches) for the 
same values of F*, /3 and 5 obtained in sections IIV Al 
and IIVBI To be more specific, the data from all the 
samples are plotted using the same couple of values for 
/3 and 6, chosen equal to the average values obtained 
in previous sections ; F* on the other hand varies from 
sample to sample as shown in the inset of Fig. [5] For a 
given disorder strength vq and Tq are constant, whatever 
the system size {N^ = 270 to 1080) and the positions of 
the pins. When the disorder strength is changed, vq and 
To change {vq = 3.3 10^"^ and To = 8 10^** for ap/a^ = 
0.35 ; Wo = 1-3 10"^ and To = 1.4 10"^ for ap/a^ = 1.05) 
but S± remains unchanged. 

The collapsing of all the available data indicates that: 
(i) the size effects are not relevant, i.e. the system studied 
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is large enough for a true critical regime to be observed 
(ii) neither the critical exponents /3 and i5 nor the scaling 
function S± depend on the disorder (strength and posi- 
tions of the pins) . Our results therefore show some degree 
of universality within our model. The question or larger 
universality classes for plastic depinning of periodic ob- 
jects with N = 2 displacement fields will be addressed in 
the discussion. 



V. DISCUSSION 

A large variety for the values of the /? exponent can 
be found in the literature, and the question of the exis- 
tence of universality classes is still open for the plastic 
depinning transition. In some cases no scaling fit can 
even been found or hysteresis suggesting first order de- 
pinning transition has been reported. Our present results 
strongly suggest a second order depinning transition with 
well identified critical exponents /3 and 6 and scaling re- 
lations exhibiting some universality with regard to the 
disorder strengths, and disorder realizations arc found. 
We therefore try to compare our results and in particu- 
lar the depinning critical exponents to those reported in 
similar studies. 

In the case of a displacement field of dimension iV = 1 
{e.g. Josephson junction arrays, XY model, metallic 
dots, CDW) in an embedding space of dimension d ~ 2, 
the value (3 — 5/3 has been predicted theoretically^^. 
Moreover, many studies are available, both experimen- 



tal and numerical, and a large variety of /? values can 
be found {e.g. 1.4 < (3 < 2.25 for arrays of metallic 
dots^, and 1.3 < /3 < 2.6 for Josephson junction and XY 
model^). 

In the case of a displacement field of dimension 
N — 2 {e.g. superconductor vortices, colloids, Wigner 
crystals), there are to the best of our knowlegde only 
a few numerical studies proposing critical exponents 
for the plastic depinning transition (see Table |I|.Ref. 
uses a iV = 2 model to describe arrays of metallic 
dots, and the nature of the disorder's randomness is 
different from the other studies (regular array of pins 
with random strength vs. identical pins with random 
positions). The same study also considers the case of 
random positions by removing random sites on the 
array, however no unique value of /? is measurable in 
that case. In Wigner crystals and colloids, (3 is found 
to be independent of the pinning strengthii^. In 3D 
superconductors {d — 3) with well defined vortex lines, 
it is shown that one can measure two critical exponents 
/? and (5, both independent of the disorder strength, and 
that a scaling relation involving these two exponents 
exists^. However, an hysteretic discontinuous transition 
is found if the superconductor is layered and the layers 
are allowed to decouple, i.e. when there are no more 
vortex lines^^. Finally in experiments on depinning in 
superconductors, /? ^ 1.2 or 1.4 is found above the peak 
effect^, where the depinning is considered to be plastic. 



Study 


d 


particle-particle 
interaction 


pinning 






Ref.2 


2 


Coulomb 


Coulomb 
random positions 


1.61 ±0.10 
1.71 ±0.10 




Ref.^ 


2 


screened Coulomb 


parabolic 
random positions 


1.94 ±0.03 




Ref.^ 


2 


Coulomb 


parabolic 
random strength 


1.94 ±0.15 




Ref. 24 


3 


modified Bessel 


gaussian 
random positions 


0.754 ±0.010 


0.754 ±0.010 


Our results 


2 


modified Bessel 


gaussian 
random positions 


1.3 ±0.1 


0.75 ±0.1 



Table I. Depinning exponents (simulation results) in the plastic case for periodic systems with a displacement field of dimension 
N = 2. dm the dimension of the embedding space. 



In the present paper we propose (to the best of our tic depinning in the case N ~ 2 and c? = 2 by study- 
knowledge) the first determination of both exponents (3 ing effective 2D vortex lines. First, we note that the 
and 5, along with evidence of a scaling relation, for plas- value (3 = 1.3 ± 0.1 is in agreement with vortex experi- 
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ments above the peak effect^. Second, we find some de- 
gree of universality witli regard to tlie pinning strength 
and different realizations of the disorder, in agreement 
with previous resultsii^. However, the variety of the val- 
ues of the (3 exponent reported above suggests that the 
type of particle-particle interaction as well as the type of 
disorder plays an important role. In the light of these 
results, it seems therefore difficult to define large univer- 
sality classes for the plastic depinning of periodic objects 
in 2D. 



VI. CONCLUSION 

In this paper, we studied the plastic depinning of vor- 
tex lattices in 2D with strong disorder. Above the pinned 
region, we find four dynamical regimes: periodic, disor- 
dered chaotic, smectic chaotic and decoupled channels. 
The first one is controlled by the finite size of the simu- 
lation box, leading to the so-called single particle regime. 
A continuous (second order) phase transition is found at 



the depinning threshold. The critical regime associated 
to the depinning transition is found to be chaotic. A 
careful analysis at T = and finite temperature allows 
an accurate determination of the critical force. Scaling 
laws for the force and temperature dependence of the 
velocity are found and two critical exponents are deter- 
mined. Both critical exponents and the scaling function 
arc independent of the disorder (strength and positions 
of the pins) in the range of parameters we have stud- 
ied, indicating some degree of universality. However, the 
comparision with similar systems may suggest that large 
universality classes do not exist for the plastic depinning 
in 2D. 
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